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ABSTRACT 

Aims. The binary system rj Carinae has completed its first 5.54 y orbit since the beginning of science operation of the 
Fermi Large Area Telescope (LAT). We are now able to investigate the high-energy 7 -ray source at the position of 77 
Carinae over its full orbital period. By this, we can address and confirm earlier predictions for temporal and spectral 
variability. 

Methods. Newer versions of the LAT datasets, instrument response functions and background models allow for a more 
accurate analysis. Therefore it is important to re-evaluate the previously analyzed time period along with the new data 
to further constrain location, spectral shape, and flux time history of the 7 -ray source. 

Results. We confirm earlier predictions of increasing flux values above 10 GeV toward the next periastron passage. 
For the most recent part of the data sample, flux values as high as those before the first periastron passage in 2008 
are recorded. A comparison of spectral energy distributions around periastron and apastron passages reveals strong 
variation in the high-energy band. This is due to a second spectral component that is present only around periastron. 
Conclusions. Improved spatial consistency with the 7 -ray source at the position of r] Carinae along with the confirmation 
of temporal variability above 10 GeV in conjunction with the orbital period strengthens the argument for unambiguous 
source identification. Spectral variability provides additional constraints for future modeling of the particle acceleration 
and 7 -ray emission in colliding-wind binary systems. 
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1. Introduction 


Five years have passed since the first reports of the unex¬ 
pected presence of a strong high-energy 7-ray source at the 


position of the enigmatic bin ary system rj Carinae (Tavani 
et al.|2009| |Abdo et al.|2010|). This binary system is located 
2.3 kpc away from Earth and comprises a massive LBV star 
and an O- or B-type companion star in an highly eccen¬ 
tric orbit (e ~ 0.9). Since the first reports, several stud¬ 
ies have tried to further investigate the system’s nonther- 
mal emission. Subsequent analyses with t he Fermi-hAT 
(Farmer et al. 2011 Reitberger et al. 2012) indicated flux 
variability in concordance with the orbital period of the 
binary system. Although the source was detected up to en¬ 
ergies of 300 GeV with Fermi-LAT, ground-based Imaging 
Atmospheric Ch erenkov Telescopes (IACT ) did not yet re¬ 
port a detection ( Abramowski et al.|2012[ ). Below the MeV 
range, e vidence remains tor a nontherrnal hard X -ray com¬ 
ponent ( Leyder et al.||2008 Sekiguchi et al. 2009), which is 
presumably linked to the emission as seen with Fermi-LAT. 
Apart from alternative suggestions (e.g., |Ohm et al. 


2010), the most common interpretation of 7-ray emission 
in The 77 Carinae binary system is a colliding-wind bi¬ 
nary scenario. In this, charged particles are accelerated 
at the shock fronts of an extensive wind-collision region 
(WCR), which forms wher e the powerful stellar winds of 
the high-mass stars collide (Eichler fc Usov|1993 |Pittard fc 


Dougherty 2006). Traveling alongside the hot, shocked gas 
in the WCR, the charged particles (leptonic and hadronic 
in nature) interact with stellar radiation fields, magnetic 
fields, and the surrounding plasma, subsequently emitting 
nontherrnal radiation via inverse Compton emission, syn¬ 
chrotron e mission, relativistic bremsstrahlung, and neutra l 
pion decay ( Reimer et al.|[2006 [Bednarek fc Pabich 


2011 ). 


Because the properties of the WCR change drastically with 
stellar separation, the ensuing nontherrnal emission is ex¬ 
pected to show strong modulation on orbital timescales for 
systems with high eccentricity such as rj Carinae. 

With the conclusion of rj Carinae’s first full orbit in 
Fermi-hAT 7-ray data and improved analysis performance, 
it is now time to perform a dedicated data analysis to de¬ 
termine whether previous indications and predictions for 
flux variability can be confirmed over one orbital period in 
7-rays. 


2. Observations with the Large Area Telescope 

2.1. Dataset 

The time range of the Fermi-LAT dataset (P7 reprocessed 
SOURCE class) that is used in this analysis spans 2024 
days from 2008 August 4 to 2014 February 18 (Mission 
Elapsed Time: 239557417 s to 414431017 s). During this 
period, the position of 77 Carinae has been monitored for a 
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total of 399 days (at energies 0.2 to 10 GeV) and 392 days 
(at energies 10 to 300 GeV). Different exposures for differ¬ 
ent energy bands are due to the energy-dependent instru¬ 
ment response functions (IRFs). To account for irregular 
exposure due to targeted observations and other deviations 
of the prevalent regular sky-coverage of the FerTTiJ-LAT, 
we considered bins of equal exposure time rather than real 
time in our temporal analysis. The data were reduced and 
an^zed using the Fermi Science Tools v9r31 pack¬ 
age To reject atmospheric 7 rays from the Earth’s limb, 
time periods when the region was observed at zenith angle 
greater than 105° were excluded and a rocking angle cut 
of less then 52° was applied. We used all remaining pho¬ 
tons with energy E > 200 MeV within a square of 21° base 
length (aligned in directions of right ascension a and dec¬ 
lination (5) centered on the nominal location of r] Car, (a, 
^) = (161.265°, -59.685°). 


differential flux. 

Source detection significance can be described by the 
likelihood test statistic value TS = —2 ln(L„iax.o/Amax,i), 
which compares the ratio of two values that are obtained by 
a maximum-likelihood procedure. Tmax.o is the likelihood 
for a model without an additional source at a specified lo¬ 
cation (the null-hypothesis), and Lmax.i is the maximum- 
likelihood value for a model including an additional source. 
Surrounding sources and background are taken into ac¬ 
count in both cases. The same notion can also be applied to 
comparing the likelihood of two source models that are de¬ 
scribed by different spectral models for a specific source of 
interest. The corresponding TS-value can be converted into 
detection significanc e a for a point source via a ~ 'sjTS 
( Mattox et al.|[l9% ). 


2.2. Likelihood analysis 

Because rj Garinae is located at low Galactic latitude within 
the projection of the Garina-Sagittarius Arm of the Milky 
Way, the diffuse Galactic 7 -ray background significantly 
affects the detected 7 -ray flux and needs to be modeled 
with care. The size of the Fermz-LAT point-spread func¬ 
tion (PSF) makes knowing the neighboring 7 -ray sources 
and their contributions equally important. The 68 % con¬ 
tainment radius corresponds to more than 2° at 200 MeV 
(on-axis). 

We performed maximum-likelihood analyses using the 
P7REP_SOURCE_V15 IRFs in conjunction with the 
Galactic diffuse model template-4years-P7-V 15-repro-v3.fits 
and the isotropic background component iso_source-v05. txt. 
The applied source model for all subsequent investigations 
includes 199 point sources that are located within a radius 
of 20 ° around the source of intere st and are part o f the 
Fermz-LAT 4-year (3FGL) catalog ( Acero et al.||2015 ). For 
each source we used the same spectral model as in the cat¬ 
alogue. The spectral parameters of all the sources within a 
radius of 10 ° were allowed to vary in the likelihood analy¬ 
sis. Most sources were modeled either by a simple power-law 
(PL) 


dN_ / E 

Ile ~ 


( 1 ) 
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Exceptions are 7 -ray pulsars, which were modeled by a cut¬ 
off power-law (GPL) function expression 


^=B ( — 
dE \Eo 


exp 


E 


( 3 ) 


Eb is the cutoff energy. In all three cases, ^ gives the 


where A is the normalization parameter and a is the spec¬ 
tral index, or by a LogParabola function, which provides 
additional information on the spectral curvature of some 
sources by introducing a break energy E^ and the parame¬ 
ters a and /? 


where B is the normalization, /3 is the spectral index and 


^ See the Fermi Science Support Centre (FSSC) website for 
details: http: //f ermi . gsf c .nasa.gov/ssc/ 


3. Analysis results 

3.1. Spatial analysis 

In the two celestial coordinate grids in Fig. the TS- 
values indicate the detection significance for a trial addi¬ 
tional point source (a PL with two free parameters) that 
was considered sequentially at each of the grid positions. 
The left figure covers the region around 3FGL J1045.15941 
(the name of the 7 -ray source that is no w associated with 77 
Garinae according to the 3FGL catalog ( Acero et al.|2015| ) 
at energies from 0.2 to 10 GeV, the right figure at ener- 
gies from 10 to 300 GeV, each representing 2024 days of 
data. The energy band was divided according to the pre¬ 
vious indication of two different emission components. The 
two TS-maps were obtained with the tool gttsmap with 
r] Carinae removed from the source model. In addition, con¬ 
fidence contours mark the 68.3%, 95.4%, and 99.7% uncer¬ 
tainty regions for the location of the source defined by the 
maximum TS value in the grids. In both energy bands the 
nominal position of rj Carinae is found within the 95.4% 
error region. For the low band it is at the border of the the 
68.3% error region. 

In the low-energy band the 7 -ray signal is centered 
on {a, (5) = (161.265°, —59.698°) with a fairly circular 
95.4% confidence region of radius r = 0.014°. In the high- 
energy band the signal is centered on a {a, 6) = (161.278°, 
—59.698°) with with a fairly circular 95.4% confidence 
region of radius r = 0.021°. Statistical limitations toward 
higher energies yield a larger confidence region for the hard 
band despite the larger PSF at lower energies. We find 
that localization uncertainties (in terms of error contour 
radii) have decreased with respect to an earlier analysis 
(Reitberger et al. 2012) by a factor of 1.6 for the low band 
and 2.3 for the high band. At the same time, the angular 
separation of the best-fit positions for the energy bands has 
decreased from 0.033° to merely 0.007°. 

The difference to the indicated catalog position of 
3FGL J1045.15941 is understood by acknowledging the dif¬ 
ference in exposure and the decomposition into low- and 
high-energy bands. 

Comparing these resul ts to an earlier localizati on study 
with 35 months of data (Reitberger et al. 2012), spatial 
agreement in both energy bands is considerably improved 
(reduced error contours). The reason is not only the better 
statistics, but also the continuously improving IRFs and 
knowledge about residual backgrounds. The spatial agree- 
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Right ascension Right ascension 

Figure 1. Test statistics maps of the celestial region around 3FGL J1045.15941 as obtained by the likelihood analysis. 
The TS-value corresponds to the likelihood of the source being located in a specific grid point. Left plot: 0.2 to 10 
GeV, right plot: 10 to 300 GeV. The confidence contours (green for the low band, white for the high band) mark the 
68.3%, 95.4%, and 99.7% uncertainty regions for the location of the source of interest. The dashed yellow circle marks 
the cataloged 95% source location region of 3FGL J1045.15941. The nominal position of ry Carinae is indicated by the 
yellow cross. 


ment between the detected y-ray emission and location of ry 
Car, as well as the agreement of the y-ray source positions 
for the consecutive energy bands, is better than it has been 
for fewer datasets. 

Performing a maximum-likelihood fit using the tool gt- 
like (dataset and source model as described above) with the 
source of interest fixed at the nominal position of ry Car, we 
determined the total y-ray flux values for the two respec¬ 
tive energy bands. The source of interest was modeled with 
a GPL for the low-energy band (0.2 to 10 GeV) and a PL 
in for the high-energy band 10 to 300 GeV). Systematic un¬ 
certainties were estimated using bracketing IRFs that take 
into account s ystematic variations of th e Fermi-hAT effec¬ 
tive area (see Ackermann et al. 2012). At the low-energy 
band, the source has a flux of F <loGev “ (1-00 ± O.OSgtat. 
^013 syst.) X 10“^ cm“^s“^ and a TS value of 3660 (~ 
60 ct). 

At the high-energy band, the source has a flux of 

LOGeV _ fa no _L. +a71 -1-0.43 \ „ in-l 


^<300GeV “ (6.03 ± _Q 07 stat. _o 37 syst.) X 10 
and a TS value of 286 (~ 17cr). 


cm ^s ^ 


3.2. Spectral analysis 

The energy spectrum of rj Carinae was determined by ap¬ 
plying the same maximum-likelihood fitting technique as 
above to individual energy bins that were chosen to be con¬ 
sistent with earlier studies ( Reitberger et al.]|2012 ). 

Figure (left) shows^ tEe spectrum tor the to¬ 
tal dataset of 2024 days and is best described by a 
SmoothBrokenPowerLaw (SBPL) function of the form 


with the fit parameters as given in Table 

Quite different spectra are obtained for a period of 
500 days centered on the apastron and periastron passages 
(phase 0.5 and 0 or 1), see Fig.(right). The choice of these 
two time intervals is a compromise between the objective 
to obtain representative spectra for the time of lowest and 
highest stellar separation and the need for sufficiently high 
(and roughly equal) statistics to allow for comparison of 
the two. Owing to the early occurrence of the periastron 
passage shortly after the beginning of LAT observations in 
2008, the respective interval includes data from the begin¬ 
ning of the full-orbit dataset (^410 days) and the end (~90 
days). The two distinct spectra indicate that the full-orbit 
spectra cannot be viewed as representative for each phase 
of the orbit. Instead, it is the average of quite distinct or¬ 
bital states. This predominantly concerns the high-energy 
regime above ~20 GeV where the difference between the 
apastron and periastron spectra is evident. For the energy 
bin 32 to 100 GeV (second from right) the periastron value 
of the differential flux is a factor of 10 higher than the up¬ 
per limit determined at apastron. The fit to the periastron 
spectrum is a GPL-fPL function (following the approach 


Reitberger et al. (20121), which serves best to describe 


the two-component spectral shape. An alternative fit to the 
periastron spectrum by a GPL model yields an unsatisfac¬ 
tory presentation of the data {x^/dof > 2.2) and was re¬ 
jected against the GPL-fPL with more than 2cr. The F-test 
between both spectral models gives a value of 3.87, which 
rejects the null hypothesis that both models fit equally well 
at 94% confidence. 


dE \EoJ 




For the apastron spectrum the second component at 
high energies disappears along with minor less significant 
alterations at the low-energy band. The favored spectral 
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Figure 2. Spectral energy distribution of rj Carinae obtained by likelihood analysis (points). The error bars are of la 
type and only reflect the statistical uncertainties. The upper limits (represented by arrows) have been determined such 
that the difference of the logarithmic likelihood values (with and without an additional trial point source at the indicated 
flux value) correspond to la. The line depicts the best fit of the spectral model that best describing the spectra. Left: 
representing the full orbit dataset. Right: representing 500 days centered on the periastron passage (blue) and the apastron 
passage (red). 


parameter 

SBPL 

full-orbit 

CPL-tPL 

periastron 

CPL 

apastron 

bin 

At (days) 

0.2 to 10 GeV 

At (days) 

10 to 300 GeV 

A 

38 ± 50 

25.2 ± 9.8 

20.4 ± 1.2 

10 "^^ MeWm-^s-B 

209 

216 

a 

1.54 ± 0.89 

2.0 ± 0.1 

2.09 ± 0.06 

2 

211 

214 

B 

_ 

11.1 ± 4.4 

_ 

10 "^^ MeV^cm-2s"^3 

195 

192 

P 

2.39 ± 0.06 

1.35 ± 0.57 

- 

4 

197 

197 

Eb 

0.60 ± 0.56 

1.5 ± 1.0 

16.9 ± 6.1 

GeV 5 

197 

192 

e 

0.28 ± 0.17 

- 

- 

6 

198 

198 


Table 1. Best-fit parameters for functions (4), (1-1-3), and g 

(3) fitted to the full-orbit, periastron and apastron spectra. g 

The scale energy Eq was set to 1 GeV. 10 


195 

203 

201 

219 


194 

202 

200 

219 


Table 2. Time bins of equal exposure as used for the tem¬ 
poral analysis. 


model now is a simple CPL function. The respective pa¬ 
rameters are listed in Table [TJ 

We studied the SEDs of the nearby pulsar PSR 


J1048—5832 as a cross-check (see comparison in Reitberger 


et al. (2012) for the same 500-day time intervals) and did 


not find any significant variations in its spectrum. Measured 
data points overlap within the Icr error bars. 


3.3. Temporal analysis 

We now investigate the flux time history of the observed 
7 -ray emission. As in the previous section, we analyzed the 
two adjacent energy intervals 0.2 to 10 GeV and 10 to 300 
GeV. 


For each time bin we fit a CPL function (all parameters 
free) to the data. Taking the average flux as determined 
for the whole dataset (see Sect. as the hypothesis for 

a non-varying source, a reduced of 0.78 for 9 degrees 
of freedom was found. This corresponds to the observation 
in SEDs in Sect. |3.2[ where no significant variability was 
detected in the low-energy part of the spectrum. Owing to 
the larger statistics, the relative errors on the flux values 
in the low-energy band are substantially smaller than those 
in the high-energy band, which, in principle, would imply 
a higher receptivity for establishing flux variability than in 
the high-energy band. 


3.3.1. Flux studies in the low-energy band (0.2 to 10 GeV) 

In Fig. left, we show the flux time history of rj Carinae 
for the energy band 0.2 to 10 GeV as obtained by likelihood 
analysis. The dataset was divided into ten consecutive in¬ 
tervals of equal exposure. The durations of each bin in real 
time are listed in Tabled 


3.3.2. Flux studies in the high-energy band (10 to 300 GeV) 

In Fig. right, we show the corresponding flux time his¬ 
tory for the high-energy band. For each time bin a PL 
function (all parameters free) was fit to the data. Again, 
ten time bins of equal exposure (as listed in Table were 
considered. The variability is far more pronounced tnan in 
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the low-energy band. The fluxes for the first two bins— 
roughly centered on the periastron passage—are about a 
factor of 3.5 higher than the fluxes for the subsequent four 
bins that cover the orbital phases before the apastron pas¬ 
sage. After apastron, the flux shows a roughly gradual in¬ 
crease until it reaches the level of the previous periastron 
flux shortly before the next periastron passage. A y^-test 
(reduced of 2.00 for nine degrees of freedom) shows that 
a null-hypothesis of observing a constant average flux (as 
determined in Sect. is disfavored at the level of 2.5cr. 
The X-ray light-curve does not resemble sinusodial modu¬ 
lation, but has a spike plus deep abso rption trough around 
periastron ( Moffat fc Corcoran|2009 1, therefore we devised 
a different statisticaf test to quantify the indicated orbital 
modulation. Under the assumption that orbital modulation 
is marked by highest 7 -ray flux aro und periastron, as for¬ 
mulated in Reitberger et al. (20121 on th^basis of three 
fifths of the ten light-curve (Ic) bins in Fig. the complete 
Ic over the full orbital phase can now be investigated.To 
avoid introducing biases, we did not select Ic bins individ¬ 
ually that may or may not relate uniquely to a periastron 
high-flux condition. Instead, we tested a multitude of Ic bin 
combinations and subsequently corrected the obtained re¬ 
sult for trials. For this we considered all fractions of the 
Ic comprising two, three, and four consecutive Ic bins that 
include the periastron passage and that are roughly sym¬ 
metric with respect to phase zero. 

The 7 -ray flux observed over the periastron passage dif¬ 
fers by 4.5 ct (post-trial significance) from the rest of the 10 
to 300 GeV 7 -ray light curve. 

The significant variability in concordance with the or¬ 
bital period at energies 10 to 300 GeV agrees with what 
we reported fro m the spectra of Sect. and also with an 
earlier analysis (Reitberger et al. 2012 1 . There are some dif- 
ferences betweeri the light curves in Reitberger et al. ( 2012 ) 
and those presented here, however. We highlight three as¬ 
pects before we discuss this: First, this analysis was carried 
out using newer versions of the LAT datasets, instrument 
response functions, and background models, which allows 
for a more accurate analysis now. Second, we changed from 
Ic bins that are equally spaced in time to those based on 
equal exposure, now aimed for minimizing effects from ex¬ 
posure differences between individual Ic bins. Last, even 
then the differences are at low significance. Comparing the 
light curves nevertheless, there is a mild indication that in¬ 
cluding the exact time of the periastron passage tends to 
produce a mildly lower gamma-ray flux than in near peri¬ 
astron bins that do not. This hint neededs to be studied 
in more detail on short orbital phase intervals. Because the 
Fermi data do not yield significant detections on distinc¬ 
tively shorter Ic bins, it might be investigated using model 
predictions. 

4. Discussion and summary 

We have shown that flux variability at energies above 
10 GeV prevails throughout the complete first full-orbit 
dataset of a 7 -ray source at the position of ij Car. Using 
a new dataset (Pass? reprocessed), revised and updated 
IRFs (P7REP_SOURCE_V15 instead of P7SOURCE_V6), 
a new source model (based on the 3FGL instead of the 
2FGL catalog), as well as updated backg round models, we 
confirmed the previously observed (see Reitberger et al. 


In addition, we tested the previously predicted increase in 
flux at high energies that had been expected to occur after 
the last apastron passage. It has been shown that the flux 
recently attained levels that resemble those before the last 
periastron passage in late 2008. Although caused by non- 
thermal processes instead of thermal ones, the high-energy 
7 -flux light curve above 10 GeV shows some general qual¬ 
itative resemblance wit h 77 Carinae’s th ermal X-ray signal 
as measured by RXTE ( Corcoran|2005 ), related to the fact 
that both depend on conditions in the binary’s WCR. 

The full-orbit dataset allowed us to compare the spec¬ 
tra centred on the apastron and periastron passages. They 
clearly deviate in the high-energy part, which can be mod¬ 
eled by an additional hard PL component emerging around 
the periastron passage. For energies from 32 to 100 GeV, 
the observed flux at periastron is more than a factor of 10 
higher than the upper limit determined at apastron. 

The phase-locked flux changes in the high-energy emis¬ 
sion component together with improved spatial agreement 
with the location of f] Garinae strengthens the case for a 
firm identification. 

We currently have two scenarios that may account for 
the spectral variation observed in 77 Carinae in the context 


7 -ray emission ( 

Eichler & Usov 

|1993 Reimer et al. 

2006 

Pittard & Dougl 

ierty||2006p. 'i'he spectral extension emerg- 


2012 1 decrease in 7 -ray flux at energies from 10 to 300 GeV. 


ing at periastron could (as proposed by Farmer et al. 2011) 
be due to an additional (fully or partly hadronic) 7 -ray 
emission component. This scenario is also plausible in the 
light of recent 3D hydrodynamical studies that take into ac¬ 
count the dynamics of the WCR and simultaneously solve 
a transport equation for shock-accelerated particles (see 
Reitberger et al. |2014b|a| ). Although not yet applied to the 
specific parameters of the 77 Carinae binary system, these 
models show that considerable variations—even a transi¬ 
tion from lepton- to hadron-dominated 7 -ray emission— 
occur with changing stellar separation in a typical colliding- 
wind binary system. It remains to be seen what the en¬ 
hanced 3D hydrodynamic models predict for the 77 Carinae 
binary systems when they are applied to its specific stellar 
and orbital parameters. However, a comparison of model 
predictions and observational 7 -ray data for this peculiar 
binary system is extremely challenging because of the com- 
plexity and turbulen ce of the wind collision region (see e.g., 
Madura et al.|[2M3 ). 

An alternative to postulating a second emission com¬ 
ponent emerging around periastron relates to phase- 
modulated strong photon-photon absorption in the complex 
ra diation fields in the vic inity of 77 Carinae. As discussed 
in iReitberger et al.| (|20I2|), the observed spectral variation 
might be explained oy the presence of hot shocked gas in the 
WCR or by spatially extended X-ray emission components 
surrounding the system and the ensuing photon-photon ab¬ 
sorption. Different properties of an absorber at periastron 
and apastron conditions could still respond to the observed 
spectra, rendering such a scenario less likely, although not 
firmly excluding it. 

Ongoing sophisticated modeling efforts may provide 
more insights into these questions. The high-energy 7 - 
ray flux up to 300 GeV around periastron raises expecta¬ 
tions that it might be detected with the lACT instrument 
H.E.S.S 2, at least around the periastron passage. 
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Figure 3. Left: Flux time history of ij Carinae in the 0.2 to 10 GeV energy band obtained by likelihood analysis. The 
lower x-axis is in units of orbital phase of the 77 Carinae binary system, the upper in units of Modified Julian Day (MJD). 
The flux error bars are of Icr type. The gray data points on the right are the phase-shifted counterparts of the first two 
data points on the left. The dashed vertical line indicates the time of the periastron passage (blue) and the apastron 
passage (red) of the rj Carinae binary system. Right: The same for energies 10 to 300 GeV . 
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